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Spectroscopic follow-up of 70 /im sources in Spitzer Wide-area 
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ABSTRACT 

We present spectroscopic follow-up observations of 70 pim selected galaxies from the SWIRE 
XMMLSS and Lockman Hole fields. We have measured spectroscopic redshifts for 293 new 
sources down to a 70 /im flux limit of 9mJy and r < 22 mag. The redshift distribution peaks 
at z ^ 0.3 and has a high redshift tail out to z = 3.5. We perform emission line diagnostics 
for 91 sources where [OIII], H/3, [Nil], Ha and [SII] emission lines are available to determine 
their power source. We find in our sample 1 3 QSOs, 1 Seyfert II galaxy, 33 star forming galax- 
ies, 30 composite galaxies, 5 LINERs and 21 ambiguous galaxies. We fit single temperature 
dust spectral energy distributions (SEDs) to 81 70 fim sources with 160 fim photometry to 
estimate dust temperatures and masses. Assuming the dust emissivity factor (/?) as 1.5, we 
determine dust temperatures in the range ^ 20-60K and dust masses with a range of 10^-10^ 
M0. Plotting these objects in the luminosity-temperature diagram suggests that these objects 
have lower dust temperatures than local IR luminous galaxies. The Herschel Space Observa- 
tory will be crucial in understanding the nature of these sources and to accurately determining 
the shape of the Rayleigh-Jeans tail of the dust SED. We then model SEDs from optical to 
far-IR for each source using a set of galaxy and quasar templates in the optical and near-IR 
(NIR) and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst and 
AGN dust torus) in the mid-IR (MIR) to far-IR (FIR). The number of objects fit with each dust 
template are: 57 Arp 220, 127 M82, 9 cirrus, 1 AGN dust torus, 70 M82 and ciiTus, 26 M82 
and AGN dust torus and 3 Arp 220 and AGN dust torus. We determine the total IR luminosity 
(Lir) in range 10^-10^^ Lq by integrating the SED models from 8 to 1000 /im. 
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1 INTRODUCTION 

The discovery of the Cosmic I nfrared Background (CIRB; 
IPuget et alJll996l : lFixsen et alj|l998l has shown that at least half of 
the energy generated by star-formation and Active Galactic Nuclei 
(AGN) in the Universe has be e n absorbed by dust and re-radiated 
in the IR dCispert et al1l2000l : iHauser & Dwekll200lh . A signifi- 
cant fraction of the background has been resolved into a population 
of IR sources, luminous (LIRG: Lm, = 10"-10^^ L©), ultralumi- 



nous (ULIRG: Lir : 



Lq) and hyperluminous (HLIRG: 



Lir > 10 Lq) infrared (IR) galaxies whose bolometric energy 
output is dominat ed in the IR (A = 8-1000 /xm) by reprocessed 
dust emission (see ISanders & Mirabel 19961: Ichary & EIb"a3l200ll : 



iFranceschini et al] 



20011) . They were first catalogued in the local 



Universe (z ^ 0.3) by the InfraRed Astronomical Satellite (IRAS) 
and shown to only contribute ~ 30% to the CIRB which suggests 
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that the majority originates from reprocessed dust emission by high 
redshift galaxies. 

The Spitzer Space Telescope jWemer et ai]l2004) has greatly 
increased our understanding of the IR Universe at high redshifts. 
We now know that there is strong evolution in the population of 
IR galaxies with LIRGs dominating the cos mic star fo rmation rate 
out to z 1 and ULIRGs out to z ~ 2 jEe Floc'h et al. 200j; 
IPerez-Gonzalez et al.ll200j : ICaputi et alj|2007l) . This evolution is 
also seen in th e more luminous s ubmillimeter galaxies (SMGs) at 
redshifts > 2 jSmail et al.ll 19971 : IScott et al.ll2002l : Iciements et all 
l2008l : lDve et al.l l2008')."These studies have highlighted that a large 
fraction of star formation and gravitational accretion is heavily ob- 
scured by dust and therefore studying the properties and the nature 
of these galaxies is crucial to understanding galaxy evolution. 

Deep surveys conducted by Spitzer, particularly in the mid- 
IR (MIR) at 24 p,n\ where the MIPS instrument is most sensitive 
have allowed us to understand the nature of dust obscured galax- 
ies. Using spectroscopic diagnostics and/or modelling the spectral 
energy distributions (SEDs) combining optical and IR photome- 
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try has shown that star-formation is the dominant process with 
the AG N component becoming significant in the more luminous 
sourc es dOenzel & CesarskvUlOodlFarrah et alll2003l : ISoifer et alj 
l2008h . However the 24 surveys are limited as the peak of the 
SED for most IR luminous galaxies is in the range of ~ 40 — 200 
^m. At high redshifts the 24 /^m channel samples shorter wave- 
length regions contaminated by polycyclic aromatic hydrocarbon 
(PAH) emission and silicate absorption features making it difficult 
to obtain reliable estimates of the bolometric IR luminosities. 

Spitzer also observed at the longer 70 and 160 wave- 
lengths and several studies have been carried out in order to char- 
acterise the 70 /^m population. Analysis of the 160 population 
is severely hindered because of the lower resolution and sensitiv- 
ity of the MIPS instrument at this wavelength. By modelling the 
IR SEDs o f 70 iim selected g alaxies in th e 0.1 < z < 2 red- 
shift range, ISvmeonidis et all bOOTl l2008h found that these ob- 
jects require a cold (dust temperature < 20K) emission compo- 
nent in the FIR. In the follow-up study, ISvmeonidis et al. I j2009h 
(henceforth S09) suggest that the MIPS 70 /xm population may be 
the missing link between the cold z > 1 submillimeter (submm) 
population detected by SCUBA and the local IR luminous galax- 
ies. Kartaltepe et al. (2010) extend their study to higher redshifts, 
z ~ 3.5 for 1503 reliable and unconfused 70 /im selected sources 
in the Cosmic Evolution Survey (COSMOS). Their analysis shows 
that the SED shapes are similar to local objects but also find evi- 
dence for a cooler component than is observed locally. Using an- 
cillary radio and X-ray data, they find that the fraction of AGNs 
increases with Lir, and nearl y 51% of ULIRGs and all HLIRGs 
likely host a powerful AGN. S vmeonidis et alj j20IOl) also arrive at 
the same conclusion although they find 23% of ULIRGs contain an 
AGN. These figures are consistent with observations of l ocal UL- 
RIGs (see .Veilleux et al, 1995 ; Sanders & Mi rabel 1 994; , Kim et al.l 
Il998l) which show that the AGN fraction increases from ~ 4 % at 
Lir = 10^° Lq to > 50% at Lir > 10^^ Lq. 

In this paper, we examine the optical and IR properties from 
a spectroscopic follow-up of 293 70 /im selected sources from 
the XMM-LSS and Lockman Hole (LH) region of the Spitzer 
Wide-area InfraRed Extragalactic survey (SWIRE) jLonsdale et al.l 
l2003 l). The present paper is the first in a series which aims to study 
the evolution of the 70 ^m luminosity function. This paper is or- 
ganised as follows: In Sect[2]we present the sample selection cri- 
teria and in Sect[3]we describe the observations and data reduction 
method. We present the main results in Sect|4]and discuss the im- 
plications in Sect|5l Throughout this paper, we assume a ACDM 
cosmology with Ho = 70kms~^Mpc~^, JIa = 0.7, and f^M = 0.3. 
All magnitudes are in the AB system unless otherwise stated. 



2 THE DATA 

In this paper we have used data from t he latest release of the 
SWIRE photometric redshift catalogue of iRowan-Robinson et al.l 
( l2008l) . The catalogue contains photometric redshifts for over 1 
million IR sources, estimated by combining the optical and IRAC 
3.6 and 4.5 /im photometry to fit the o bserved SEDs with a 
combination of galaxy and AGN templates jBabbedge eTZIbOO^ ; 
iRowan-Robinson et al.ll2()05i) . Specifically, we make use of multi- 
wavelength data from XMM-LSS and LH-ROSAT regions of the 
SWIRE survey. 



2.1 Infrared data 

The SWIRE survey jLonsdale etakl l200l liooj) is one of the 
largest Spitzer legacy programmes covering 49 deg^ in six differ- 
ent fields with both the IRAC and MIPS instruments. Typical 5a 
sensitivities are 3.7, 5.3, 48 and 37.7 ^Jy in the IRAC 3.6, 4.5, 5.8 
and 8 ^.m bands. For MIPS the 5a limits are 230A (Jy, 20 mJy, and 
120 mJy at 24, 70 and 160 /xm jSurace et al.ll200a) . The IRAC and 
MIPS data were processed by the Spitz er Science Centre (SSC ) 
and sources extracted using SExtractor (Bertin & Amouts 199^. 
Full d etails of the SWIRE data release can be found in Surace et al] 
( l2005l) . 

The final data product consists of a bandmerged IRAC and 
MIPS 24 pm catalog and single-band catalogs at 24, 70 and 160 
pm. The IRAC and MIPS 24 pm. catalog consists of sources de- 
tected with a signal-to-noise ratio (S/N) > 5 in one or more IRAC 
bands and their 24 pm associations with a S/N > 3. The MIPS 
70 and 160 pm sources single-band catalogs were matched to the 
IRAC and MIPS 24 pm bandmerged catalogs to produce a SWIRE 
IRAC and MIPS 7-band catalog. Thus, almost all of the 70 and 160 
^m sources are also detected at 24 pm. 



2.2 Optical data 

Optical photometry is available for > 70% of the SWIRE area in at 
least three of the U, g', r\ i' and Z photometric bands. Spitzer- 
Optical cross-identifications (XID) was carried out between the 
optica l and the IRAC-24 ^m catalogs using a searc h radius of 
1.5" (Row an-Robinso n et al.'2005l; [Surace et al. 2005). The cross- 
identification process ensured that each SWIRE source only had 
one optical m atch. Comple t eness and reliability of the XID was in- 
vestigated by ISurace et al. I j2005l) . which showed that essentially 
the 5p;/zer-optical XIDs are essentially 100% complete. The re- 
quirement that appears to eliminate spurious sources effectively 
and give a high-reliability catalog i s that the source be detected 
at both 3.6 and 4. 5 pm at S/N ^ 5 dRowan-Robinson et ani2005l ; 
ISurace et al.ll2005l) . 



2.2.1 XMM-LSS 

The XMM-LSS centred on a = 02''21™20=, S = -04'*30™00' covers 
9.1 deg^. Optical photometry is available for 6.97 deg^ of XMM- 
LSS, which was observed as part of the Canda-France-Hawaii Tele- 
scope Legacy Survey (CFHTLS) in the ii*, g , r , i' , and z' bands to 
magnitude (Vega, 5a for a point like object) limits of: 24.9, 26.4, 
25.5, 24.9 and 23.4 res pectivelfl The photometry was obtained 
from lPierre eTai] j2007h . 



2.2.2 LH-ROSAT 

The LH-ROSAT field is smaller part of the larger LH region cen- 
tred on a = I0''52™43", 5 = -l-57'^28'"48' and covers 0.25 deg^ 
Optical photometry in the U, g\ r' and bands were obtained 
using the MOSAIC camera on the 4m-Mayall Telescope at Kitt 
Peak National Observatory (KPNO). The 5a limiting magnitudes 
(Vega) a re 24.1, 25.1, 24.4 , and 23.7 in the four bands for point-like 
sources dBerta et alj2007t) . 



1 See |http://www.cfht.hawaii.edu/Science/CFHTLS/| 
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2.3 Sample selection 

We selected 70 /im sources with r < 22 mag, S70 > 9mJy and 
SNR > 5 at 70 for the spectroscopic follow-up. We found 1553 
70 fim sources in XMM-LSS satisfying the selection criteria, of 
which 22 had a measured spectroscopic redshift. In LH-ROSAT 
58 sources were found, 13 of which also had measured spectro- 
scopic redshift. Thus, our final sample contained 1553 and 45 70 
/im sources in XMM-LSS and LH-ROSAT respectively. 

2.4 Fibre configuration 

Fibre configuration was carried out using the AF2_COnfigur43 pro- 
gramme. In total there are 150 science fibres and 10 fiducial fi- 
bres which are used to field acquisition and guiding. On average 
fibres were placed on ~ 70 primary targets (70 /im sources) and at 
least 4 standard stars per mask. In addition 10-20 fibres were used 
for simultaneous sky observations which are referred to as sky fi- 
bres. Spare fibres were placed on 24 fim sources in particular high 
redshift (zpi, > 3) candidates in the XMM-LSS and LH-ROSAT 
fields, the results of which will be presented in a future paper (Hyde 
et al. in prep). Note that due to fibre collision, it was not possible to 
utilise all the science fibres. 



3 OBSERVATIONS AND DATA REDUCTION 

The targets were observed using the Wide Field Fibre Optical Spec- 
trograph (WYFFOS) and the robot positioner Autofib2 at the Ob- 
servatorio del Roque de los Muchachos, using the 4.2m William 
Herschel Telescope (WHT). We carried out the observations over 
2 runs totalling 6 nights between November 2008-2009 using the 
small fibre bundles (fibre diameter 1.6 arcsec) and the R3I6R grat- 
ing covering the spectral range ~ 4000-9500A and a spectral res- 
olution of 8.1 A. Each fibre configuration was observed for 3 hours 
giving a S/N of ~ 5ct per resolution element for sources with 
r = 22 mag. Bias, and dark frames were obtained at the start and 
end of the night. Dome flats and arc lamps (helium and neon) were 
obtained at the start and end of each observed field. Standard stars 
were also observed at the start and end of the night for flux calibra- 
tion. 

The spectra were reduced and calibrated using standard IRAF 
routines and the WYFF0S/AF2 reduction tool made available on- 
line. The tool is based on the dofibre package which was used to 
remove cosmic-rays, subtract the bias and dark frames, flat-field, 
and wavelength calibrate, producing ID spectra. The sky fibres (re- 
ferred to as sky spectra) used for the sky correction procedure were 
median combined to create a master sky spectrum. A correction 
factor was then manually determined by comparing the brightest 
sky lines (the OI 5570A or the NaD 5890A atmospheric emission 
lines) in the target and sky spectra. The sky subtracted target spec- 
tra (referred to as object spectra) were determined by subtracting 
the master sky spectrum corrected for the factor from the target 
spectra. All object spectra were examined visually to validate the 
sky subtraction accuracy. Standard star frames were reduced using 
the same method and the object spectra calibrated using the IRAF 
standard and calibrate packages and atmospheric extinction cor- 
rected for the La Palma site. Redshifts were measured using the 
IRAF package rvsao and by visual inspection of the sky subtracted 
and calibrated ID spectra by identifying emission and absorption 

^ See |http://www.ing. iac.es/Astronomy/instruments/af2/configuration.html"] 



lines. A summary of the observation log and spectroscopic redshift 
success is given in Table[3] 



4 RESULTS 

We observed a total of 477 70 fim sources in the SWIRE XMM- 
LSS and LH-ROSAT fields and measured spectroscopic redshifts 
for 293 sources (280 in XMM-LSS and 13 in LH-ROSAT), giving 
a success rate of ~ 61%. The first observation run was severely 
affected by high humidity which resulted in poor quality spectra for 
the first two fields which explains the low redshift success rate (~ 
36%). For the second observing run, we have a measured redshift 
success rate of ~ 71%. Redshifts for some objects could not be 
measured due to sky subtraction residuals resulting from poor OH 
emission line subtraction in the wavelength region ~ 7000-9500A. 
A table of our spectroscopic redshift catalogue is provided as online 
material. 

For each measured redshift, we have assigned a redshift qual- 
ity flag. The flag has values between 0-3, with classification as fol- 
lows: zflag = 0, highly reliable redshift with three or more spectral 
features; zflag = I: redshift with at least two emission lines; zflag 
= 2, redshift from one emission line and consistent with the photo- 
z estimate; zflag = 3, redshift from absorption features only and 
consistent with photo-z estimate. Our sample contains 194 objects 
with zflag = 0, 21 with zflag = I, 59 with zflag = 2 and 19 with 
zflag = 3. The redshift distribution and the magnitude distribution 
of the sources in our spectroscopic sample are shown in Fig[T] The 
redshift distribution peaks at 2; ~ 0.3 with a high-redshift tail out 
to z — 3.389. We show in Fig|2]and Fig[3]a selection of optical 
spectra for a selection of 70 and 24 /im sources. 



4.1 Pliotometric redsliift comparison 

We have compared our spectroscopic redsh ifts with the photomet- 
ric redshift estimates using the ImpZ code JBabb edge et a l .1120040 
in the photometric redshift catalogue of iRowa n-Robinso n et al.l 
('2OO8). The reliability of the photometric redshifts are measured 
via the fractional error Az/(1 + 2) for each source, examining the 
mean error Az/{1 + z), the total rms scatter CTz and the rate of 
'catestrophic' outliers, 77, defined as the fraction of the full sample 
that has 1A2/(1 + z)\ > 0.1. 

We show the comparison of spectroscopic redshift 
(log (1 + 2sp)) against photometric redshift (log (1 + Zph)) 
in Fig|4] for objects which have a reliable photometric redshift, 
that is objects with reduc ed (x^h) < 10 in th e photometric 
redshift catalogue of Row an-Robinson et al.l ( |2008|) . The symbols 
represent the redshift quality flag while the QSOs (blue triangles) 
are separated as the photometric redshift solutions are expected to 
be less accurate for AGNs than galaxies. For the 70 nm sample 
(right panel Fig|4l, the total rms scatter, cr^ was 0.14, the mean 
error was 0.035, and the rate of catastrophic outliers rj was 5.7%. 
We find in our 70 /im sample 14 catastrophic outliers, 6 of which 
are QSOs and the rest have a reliable redshift with zflag = 0. The 
24 /im sample is characterised by = 0.13, the mean error = 0.42 
and Tj = 25.8%. For QSOs as stated, the photometric codes are 
expected to perform less reliably possibly because of photometric 
variability. For the other outliers the causes may be poor optical 
photometry or photometric redshift aliasing, but overall the low 
catastrophic outlier rate is encouraging for use of photometric 
redshifts in future analysis. 
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Figure 1. Left: Spectroscopic redshift distribution. Right: r magnitude distribution. 
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Figure 2. Selected example optical spectra of 70 ^m sources observed at WHT using AF2AVYFFOS. 
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Fibre configuration 


Observation Date 


a (J2000) 


5 (J2000) 
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Table 1. Summary of the observation log. The RA and DEC are the centre of each mask. Column 5 lists the number of 70 fim sources targeted and column 6 
is the number of successful spectroscopic redshift measured in each pointing. 
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Figure 3. Selected example optical spectra of 24/xm sources observed at WHT using AF2/WYFFOS. 
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Figure 4. Comparison of spectroscopic against pliotometric redshift for objects witli X^cd ^ from the catalogue of lRowan-Robinson et al.l )2008l) for 24 
(left panel) and 70 ^m (right panel) sources. The two solid lines are the 10% error boundary in log(l + z). The data are plotted according to the redshift 
quality flag, (zflag; see text in Sect|4j, where red crosses are objects with zflag = 0, magenta pluses are objects with zflag = 1, purple crosses are objects with 
zflag = 2 and green circles are objects with zflag = 3. Blue triangles are the QSOs identified from spectra. 



4.2 Emission line diagnostic 

We have classified narrow emission line galaxies and determined 
the dominant power s ources using emission line diagnostics or BPT 
( lBaldwinetij]|l98ll) diagrams. The diagrams are based on opti- 
cal emission line ratios such as [01Il]/H/3, [NII]/Hq, [SII]/Ha and 
[OII]/Hq:. We have used a combination of [OIII]/H^, [NII]/Ha and 
[SII]/Ha emission line ratios as these ratios are insensitive tored- 
dening because they are close in wavelength. iKewlev et alj bOOlh 
derive a new theoretical 'maximum starburst line' on the BPT di- 
agram to separate AGNs and star forming galax ies, with galaxies 
lying a bove the line likely to be AGN do minated. iKauffmann et aU 
( l2003h modified the lKewlev et al.l ( 1200 ll) scheme by including an- 
other line which includes objects whose spectra contain both AGN 
and star formation. Here we use the classification scheme presented 
in|Kewley et al. ( 2006) to classify objects in our spectroscopic sam- 
ple. 

To produce the BPT diagrams we have measured fluxes for the 
H/3, [OIII], [Nil], Ha and [SII] emission lines for 91 objects which 
have good quality spectra. We show the redshift distribution for 
these sources in Fig|6l which are mostly low redshifts (z < 0.4). 
The emission line fluxes were measured using the splot routine 
in IRAF by interactively fitting a Gaussian function to each line. 
The BPT diagrams, [NII]/Ha versus [OIII]/H^ and [SII]/Ha ver- 
sus [0III]/H/3 are shown in FigH] The blu e line is the pure star 
formation curve of IKauffmann et ah 1 2003b and the red lines are 
the extreme starburst lines of lxewlev et alj ( l200ll) . Objects which 
lie below the blue curve in the [NII]/Ha versus [0III]/H/3 and the 
red curve in the [SII]/Hq versus [OIII]/H/? are star forming or HII- 
region galaxies. Composite galaxies lie in the AGN-HII mixing re- 
gion between the blue and the red curve in the [NII]/Ha versus 
[OIII]/H/3 diagram. The black line in the [SII]/Ha ver sus [OIII]/H/3 
diagra m is the Seyfert-LINER separation curve of iKewlev et al.l 
( |200lh . Objects which lie above the black and red curves are clas- 
sified as Seyfert and objects that lie above the red curve and below 
the black line are LINERs. Ambiguous galaxies are those which lie 
above the extreme star forming line in one diagram and lie below 
in the other diagram. 




Figure 6. Spectroscopic redshift distribution for 9 1 objects used in to per- 
form emission line diagnostic. 



Using the emission line diagnostic we found in our sample: 
1 Seyfert II galaxy, 34 star forming galaxies, 30 composite galax- 
ies, 5 LINERs and 21 ambiguous galaxies. In addition 13 QSOs or 
Seyfert I (Type I AGN) were identified from their broad line spec- 
tra. We find a large number of ambiguous sources which is most 
likely due to weak [SII] emission line. 



4.3 Colour diagrams 

We use the S WIRE data with t he AF2AVYFFOS spectroscopy to 
reproduce the lLacv et al.l ( l2004h IRAC colour-colour diagram. This 
diagnostic requires the sources to have detections in all four IRAC 
bands, 3.6, 4.5, 5.8, and 8 ^m. We show in Fig|7]the IRAC colour- 
colour plot for 233 sources in our sample with detections in all the 
IRAC bands. Objects that lie within th e black wedge are expected 
to be AGN dominated as suggested bv lLacvetid]j2004h . We also 
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Figure 5. Lef t: [NII]/Ha versus [OIII]/H/3 diagnostic diagram. The blue line is the lKauffmann et al pure star formation line an d the blue lin e is the 

iKewlev et alj j200lll exti'eme stai'burst line. Right: [Sn]/Ha versus [OIII]/H/3 diagnostic diagram. Red line is the extreme starburst of Kewlev et alj fcoOll) 
is used to separate the star forming Seyfert galaxies. The black solid line is the Seyfert/LINER separation from Kewley et al. 12006). In both the panels the 
red triangles are pure star forming galaxies, blue circles are c omposite galaxie s, green squares are Seyferts, orange stars are LINERs and purple asterisks are 
objects classified as ambiguous. The dashed lines represent the lHo et aljil997l) classification scheme. A typical eiTor bar is shown on the top right hand corner 
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Figure 7. IRAC colour-colour plot iLacv et al.l2004) for 233 70 /im sources 
with detection in all IRAC bands. Cyan diamonds are QSOs, red triangles 
are star forming galaxies, blue circles are composite galaxies, green squares 
are Seyferts, orange stars are LINERs, pink asterisks are ambiguous galax- 
ies and grey plus signs are galaxies without a spec troscopic classifica tion. 
The black wedge is the AGN region as defined bv lLacv etal] i2004 . We 
also plot the co l our tra cks between z = — 4 using the template SEDs of 
ICharv&Elb"a3)200lh for different types of IR galaxies. The black circles 
represent 2 = 0. Typical error bar is shown at the bottom light corner. 



pl ot the colour tracks be tween z = — 4 using the template SEDs 
of lCharv&Eiba3 ( l200ll) . 

As expected we find that most of the broad-line objects (10/12) 
lie within the AGN wedge but we also find that 2 of the broad line 
AGNs occupy the same parameter space as the star forming galax- 



ies. This is consistent with their IR SEDs which show that the star- 
burst component is the dominant process in the IR. The star form- 
ing galaxies occupy the top-left above the AGN wedge which is 
in agreemen t with previous work jDonlev et al.l200^ : lTrichas et al.l 
l2009Ll20IOh . The narrow-line AGN lie within the AGN wedge. The 
majority of the LINERs, composite and ambiguous galaxies oc- 
cupy the same region as the star forming galaxies which implies 
that these objects are most likely starburst dominated in the NIR. 

In the left panel of Fig[8] we plot the ratios of MIPS 24 and 
70 /im flux densities as a function of redshift for all sources in our 
spectroscopic sample and colo ur tracks using the SED templates of 
iRowan-Robinson et alj j2004h (The SED models are described in 
more detail in Sect |4.4| >. The S70/S24 ratio can be used to distin- 
guish between star forming galaxies and AGN dominated sources 
except for sources with strong PAH emission features redshifted 
into the 24 fim band (e.g. M82 at 2 ~ 2) a nd for AGNs with 
strong silicate absorpt ion features at z ~ I. 5 jpraver et alj|2006l : 
[elements et alj2010al) . In this diagnostic, the star forming galaxies 
are expected to have high S70/S24 ratio, whereas AGN dominated 
sources are expected to have low S70/S24 ratio. The majority of the 
sources in our sample have observed ratios within the range of val- 
ues between the Arp220 starburst and cirrus tracks which indicates 
that a mixture of starburst and cirrus or a pure starburst or cirrus 
component dominates the FIR emission. 

More recently, iMuUanev etafl feoioh analyse FIR properties 
of X-ray sources using ultradeep 70 and 24 /im Spitzer observations 
to show that the 70/24 fim flux ratio can discriminate an AGN or 
starburst dominated IR SED for a given AGN. In our sample only 4 
of the 13 QSOs appear to have flux ratios similar to AGN dust torus 
with rest of the QSOs having ratios similar to the starburst com- 
ponent, suggesting that the IR SED is starburst dominated. SED 
models (see Sect |4.4t of our QSOs shows that most of them have a 
starburst dominated IR SED. 

We show in the right panel of FigjS] the observed 70 and 160 
Hm flux ratios for 83 sources with 160 ^m detection, which are 
a rough estimate of the dust temperature for the FIR peak emis- 
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Figure 8. Left: Observed 24 and 70 /im flux ratios as a function of redshift for all 70 /tm sources. Cyan diamonds ai'e QSOs, red triangles are star forming 
galaxies, blue circles are composit e galaxies, green square s are Seyferts, orange triangles are LINERs and grey plus signs ai'e galaxies without a spectroscopic 
classification. Colour tracks usine lRowan-Robinson et aU (2004) SED templates. Right: Observed 70 and 160 /.tm flux ratios as a function of redshift for 83 
70 [im sources. Purple diamonds are QSOs, red triangles are star forming galaxies, blue circles are composi te galaxies, green squares are Seyferts, orange 
triangles aie LINERs and grey plus signs are galaxies without a spectroscopic classification. Colour tracks usine lRowan-Robinson et al] )2004h SED templates. 
Note that the AGN dust torus colour track in the right panel lies below the origin and therefore it is not shown. Typical eiTor bars are shown at the bottom right. 



sion. The S160/S70 ratio remains roughly constant across the red- 
shift range an d is consistent with be ing cirrus or starburst domi- 
nated. S09 and lKartaltepe et aP fcOlQI compare their 70 /im sam- 
ple with local sources in the S160/S70 — z plot and find a large 
fraction of galaxioes that are colder than their local equivalent. In- 
deed our analysis in Sect |4.5.f| shows that these sources do have 
dust temperatures that are colder than local IR galaxy population. 



4.4 Spectral energy distribution 

We mode l the SEDs for our 294 sources follo wing the method de- 
scribed in lRowan-Robinson et alJ bOOSLbOOSh used for the SWIRE 
photometric redshift catalogue. We have for each source photom- 
etry in at least 3 of the 5 optical U, g, r, i and Z bands and IR 
photometry from Spitzer IRAC 3.6 - 8 fim and MIPS 24 - 160 fim 
bands. The SED fitting follows a two-stage approach, by first fit- 
ting the optical to near-IR (U to 4.5 /im) SED using the six galaxy 
and three AGN templates used by Rowan-Robinson et al. ( 20q3). 
We then calculate the IR excess by subtracting the galaxy model 
fit from the 4.5 to 24 fim data. We then fit the IR excess, 70 and 
16 fim (for 83 70 fj.m sources) data points w ith the IR template 
of iRowan-Robinson et al] ( |2004 l20oi |2008|) . The IR templates 
are derived from radiative transfer models dependent on interstel- 
lar dust grains, the geometry and the density distribution of dust. 
The IR templates are: 1.) IR 'cirrus': optically thin emission from 
interstellar dust illuminated by the interstellar radiation field; 2.) an 
M82 starburst; 3.) a more extreme Arp220-like starburst and 4.) an 
AGN dust torus. We also allow the sources to be fit by a mixture of: 
1.) M82 starburst and cirrus, 2.) M82 starburst and AGN dust torus 
and 3 .) Arp220 and AGN dust torus to properly represent the IR ex- 
cess ( IRowan-Robinson & Crawford 1 19891 ; iRowan-Robinson et al.l 
l2005h . 



We show example SED fits for some of our sources in Fig|9] 
The IR SED fits suggest that 77 sources require the cold cirrus 
component, with only 7 fit with a pure cirrus template while 286 
sources in our sample require a starburst component. 32 sources 
require an AGN dust torus component, with only 1 fit with a pure 
AGN dust torus template, 28 of which occupy the AGN parame- 
ter space in Fig[7] including 10 spectroscopically identified AGNs 
(9 QSOs and 1 Seyfert). The SED fits imply that energetically, 
the 70 fim sources are powered by star formation and the AGN 
dust torus makes very little contribution to the IR emission. For 
sources fit with an AGN dust torus and an M82 starburst mixture, 
the MIR emission is AGN dominated with the FIR emission pow- 
ered by the starburst component. This finding is to be expected as 
our selection at 70 ftm directly implies that the FIR emission is 
dominated by star formation as opposed to AGN which peak at 
much shorter wavelengths ( NIR/MIR) and are weak FIR emitters 
( lAlonso-Herrero et ali2006l) . 

For the 14 AGNs (13 QSOs and 1 Seyfert) in our sample we 
find 7 (6 QSOs and 1 Seyfert) are fit by a mixture of AGN dust torus 
and M82 starburst, 2 QSOs fit with a mixture of Arp220 starburst 
and AGN dust torus, 1 QSO fit by pure AGN dust torus, 2 QSOs 
fit with M82 starburst template and 2 QSOs are fit with an Arp220 
starburst SED i.e. no AGN dust torus. The AGN dust torus com- 
ponent in these sources typically contributes --^ 40 - 100% to the 
total IR emissions. All galaxies classified as LINERs are fit with 
the M82 starburst template. Only 1 of the composite galaxies re- 
quires the AGN dust torus component, while the rest are fit with a 
pure starburst or a mixture of M82 and cirrus. Of the 21 ambiguous 
galaxies, 2 are fit with a mixture of AGN dust torus and M82 star- 
burst whilst the rest are fit with a pure starburst or M82 starburst 
and cirrus mixture. 3 of the galaxies identified as star forming from 
the BPT diagrams require an AGN dust torus component, while 
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Figure 9. SED for a sample of 8 sources at a variety of luminosities. The black lines are the best-fit template SEDs, dashed green line is the M82 starburst, 
dot-dash blue line is cirrus and orange dot-dash line is AGN dust torus component of the IR SED. All photometry has been redshift-coiTected to 2 = 0. The 
errors bars are smaller than the size of the red points. 



only 1 is Arp220 like starburst with the rest fit with a pure M82 
starburst or M82 and cirrus mixture. 12 of the unclassified sources 
require the AGN dust torus component to model the IR emission. 

We also have in our sample 83 70 /im sources that have a 160 
^m detection, with IR luminosities in the range ~ 10*'® - lO^"*'^ Lq 
and redshift range 0.021-1.369. Of these sources 49 are fit with a 
cool cirrus component and have typical luminosities of ~ 10^^ Lq. 
This shows that a large fraction of the 70 /im sources with 160 /im 
detection have significant cold dust that contributes ~ 40 - 100% to 
the total IR luminosity. Using SED fitting S09 found that the major- 
ity of the 70 ^m sources in their sample peak at longer wavelengths 
(A > 90 fim) which they interpret as evidence for evolution in the 
cold cirrus component from 2: ~ to higher redshifts. However, at 



high redshift the 160 data point probes shorter wavelength (80 
/im at 2 = 1) and therefore the peak of the emission is not well con- 
strained and a wide range of SED templates can fit the data points. 
In order to properly understand the role of cold dust we require 
observations across a wider range of FI R and submm wave lengths 
which the Herschel Space Observatory 1 Pilbratt et al.ll2O10l) is well 
placed to do. In fact Rowan-Rob inson et al.l i201(] l) have shown the 
need for cirrus templates with colder (T = 10 - 20K) dust to model 
the SEDs of Herschel sources. In Sect l4.5l we fit single tempera- 
ture modified blackbody SEDs to these sources and study in more 
detail the link between 70 /im sources and SMGs. 



We summarise the results of the SED fitting in Table |4.4l and 
show in Fig[TO]the histogram of the total IR luminosity (calculated 
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Star-forming 
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44 


86 


4 


40 


15 
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Table 2. Summary of the number of objects fit with the 7 IR SEDs separated by optical classification. 



by integrating the SED models from 8 to 1000 /im) separated by the 
best fitting IR template. Overall, ~ 96% of the sample requires a 
starburst component, with ~ 76% represented by M82 and ~ 20% 
by Arp220-like respectively. The M82 starburst template peaks at 
a shorter wavelength in the FIR and has a less prominent silicate 
absorption feature indicating less extreme obscuration towards the 
starburst region. The Arp220-like template on the other hand is 
characterised by a more prominent silicate absorption feature due to 
the higher optical depth resulting in larger obscuration. Moreover, 
the emission peaks at longer wavelengths because an additional 
cold cirrus component is needed to match the submm emissions. 
As noted above, significant fraction of 70 /xm sources with 160 
detections are fit with the cirrus template (either as a mixture or 
pure) which is needed to account for the IR emission longword of 
100 /xm. Our results suggests that in general the IR emissions in 
the 70 population is largely powered by star formation, with 
the cold cirrus component making a significant contribution. 

We show in the left panel of Fig[TT] Lm against redshift colour 
coded by the best fitting IR template. Our sample consists of 7 nor- 
mal galaxies defined as Lir < 10^" L©, 93 starburst galaxies (Lir 
= 10" - 10" Lq), 151 LIRGs, 33 ULIRGs and 10 HLIRGs. We 
find that the fraction of sources requiring an AGN dust torus com- 
ponent in the SED models increases with the total IR luminosity 
from 1.1% at Lm < lO" to 10% for LIRGs, 27.3% at ULIRGs 
and 70% at HLIRGs. This trend is also seen with the spectroscop- 
ically identified AGNs where we find 2 LIRGs, 5 ULIRGs and 7 
HLIR Gs. T his finding is consist ent with the study of IVeilleux et al.l 
l ll995h and lVeilleuxetd] ( IT999i) , whose study of IRAS selected 
galaxies shows that the A GN fraction increases w ith Lir . Our re- 
sults are in agreement with'Kartaltepe et al.'feOKf) who find in their 
sample of 70 /im sources that the AGN fraction increases from 2% 
for Lir < 10" Lq, to 10% for LIRGs, 51% fo r ULIRGs and 97% 
for HLIRGs. Similarlv lSvmeonidis et al.l ( 1201 Oh find in their 70 ^m 
sample, AGN fractions of 0% for Lir < lO" Lq, 1 1% for LIRGs 
and 23% for ULIRGs, being consistent with our results. 

In the right panel of Fig[TT] we show the distribution of total 
IR luminosity as a function of redshift coded by optical spectral 
type determined in Sect |4.2| Star forming and composite galaxies 
span a range of IR luminosities, ~ 10^ - 10^^ Lq, with most clas- 
sified as starburst (Lir = 10^*^ - 10^^ L©) or LIRGs. Ambiguous 
galaxies also span a wide range of IR luminosities, ~ 10*'^ - 10^^ 
Lq . Most of the LINERs are seen near the starburst and LIRG clas- 
sification line. The only Seyfert II object is a LIRG while Type I 
AGNs are the most luminous sources in the sample and are typi- 
cally ULIRGs and HLIRGs. Our results are consistent with previ- 
ous studies, which suggests AGNs are more likely to be found in 
the most luminous of the IR galaxy population. We note that the 
source at z = 3.374, which has an extreme IR luminosity (Lir = 
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Figure 10. Histogram of Lir separated by the best fit IR template. Red = 
m82 starburst only, cyan = M82 and cirrus mixture and black = AGN dust 
torus only 100, green = cirrus only, blue = Aip220 starburst only, orange = 
M82 and AGN dust torus mixture and purple Arp220 and AGN dust torus 
mixture. 



14.6 Lq) is consistent with the source being a QSO. However, the 
derived IR luminosity may be unreliable as the FIR SED of this 
source is not well constrained since at z = 3.374, the 70 ^m band 
corresponds to ~ 15 /xm in the rest frame. We will be able to obtain 
a better estimate of Lir for this and other sources by constructing 
full SEDs from the optical to submm using observations from the 
Herschel Multi-tiered Extragalactic Survey (HerMES) programme 
(Oliver et al in prep.). 



4.5 Dust properties 

In this section we select the 70 ^m sources with 160 ^m detection 
at z < 1.2 (to ensure we are still sampling the peak of the FIR SED, 
which at = 1.2 160 ^m corresponds to~ 73 ^tm) to model the 
FIR SED for 81 objects by a standard single temperature modified 
blackbody given by: 



(1) 



where Td is the dust temperature and j3 is the dust emissivity, cor- 
responding to the slope of the Rayleigh- Jeans tail. This model is a 
simple approach as the parameters that define the SED do not ac- 
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Figure 11. Left: Total IR luminosity (Lir) as a function of redshift for all 293 sources in SWIRE XMM-LSS and LH-ROSAT coded by best fit IR SED type. 
Right: Total IR luminosity (Ljr) as a function of redshift coded by spectral type. Typical eiTor bar of 0.1 dex is shown at the bottom right. 



count for the detailed geometrical mix of dust grains at different 
temperatures in the interstella r medium (ISM). Nevertheless sev- 
eral s t udies (lHildebrandll983l:lDunne et a l."20QO': 'Punn e & Ealesl 



200 ll; iKlaas et al 



2001 



Bendo et alj"|2 003; Blain et aiT T2003l : 



Vlahakis et al .I2OO5I " [elements et al .l2010bH have used this method 



to fit the FIR to submm SED. 

We use the same method used for the stu dy of SCUBA Lo- 
cal Univer se Galaxy Survey (SLUGS) s ources (*Dunne et al."200tf; 
IVlahakis e t al. 2005) and local ULIRGs jClemen ts et al. 2010b) by 
calculating the dust temperature using a minimisation method 
and fixing the emissivity parameter, /3 to 1.5 1 Chapman et al.l200^ : 
iFraver et al.l2006l : lKovacs et alj2006l : ICoppin et al.l2008l) . We also 
calculate the dust masses in these objects using: 



Mi = 



ni{u)B{u,Td) 



(2) 



where Sieo is the observed 160 /im flux, D_l is the luminosity dis- 
tance to the object and Hd is the dust opacity coefficient. T he value 
of Kd is 1.139 m^ kg^ at 160 /im is interpolated from iDraind 
l l2003h . Note that our derived dust masses are crude estimates be- 
cause of the uncertainties in K.d{i^) and the dust temperatures, re- 
sulting from the complexity of the properties of interstellar dust 
(see review by Draine 2003). 

Our sample of 82 sources is characterised with mean and stan- 
dard deviation for the best fitting dust temperatures, T^ = 32.7 ± 
6.7K and dust mass range, Md ~ 10^7 - 10*. 8 Mq for /3 = 1.5. 
Our dust temperature and dust mass distributio n shown in Fig|12l 
are consistent with previous studies, [praver et al. (2006) estimate 
dust temperatures, T^ ~ 30 ± 5 K for 70 and 160 /im detected 
sources fro m the Spitzer extra galactic First Look Survey (xFLS) 
for (3 = 1.5. lDunne etal.l ( I2OOO ) SLUGS sample is characterised 
by dust temperatures, Td = 35.6 ± 4.9 and dust mass range, Md = 
10^' - 10^ Mq. Dunne et al. (2000) also fit the emissivity parameter 
finding /3 = 1.3 ± 0.2. If we adopt /3 = 1.3, we determine slightly 
warmer dust temperatures, Td = 33.9 ± 7K. 

In our analysis w e have assumed the dust emissivity to be 
1.5 but other studies dDunne & Ealesl I2OOII : iKlaas et al] I2OOII : 



I Vlahakis et ai]|2005l) have shown that the FIR-submm SED may 
be better represented by two dust component fits. This model con- 
sists of dust at two different temperatures, a cold (T^ ^ 25K) and 
a warm component (Td > 30K) and the emissivity parameter fixed 
at 2. However these studies have been carried out on local sources 
for which they have observations across a wide spectral range from 
FIR (IRAS, ISO and Spitzer photometry) to submm (SCUBA pho- 
tometry). For our sources however, the MIPS 160 ^m provides the 
longest data point and therefore it's not possible to fit a multicom- 
ponent dust SED to our sources. Even constraining the dust emis- 
sivity parameter, which varies between 1-2, is not possible as this 
requires data points across a range of wavelengths from the FIR to 
submm. 

By fixing j3 = 2, our sample is characterised by Td = 30 ± 6 
K, so increasing the emiss ivity results in slightly lowe r dust tem- 
peratures. Recent study by IPlanck Collaboration! l l20ni) of Planck 
sources has been crucial in resolving the degeneracy between jS and 
Td for nearby galaxies. Their results confirm cold dust is seen in 
Planck sources, which have median Td = 26.3 K with temp eratures 
ranging from 15 - 50 K and /3 = 1.2. lAmblard et alj ilQltX) find for 
Herschel-ATLAS sources P = 1.4 ± 0.1 which is consistent with /3 
= 1.3 for SLUGS sources founne et al]|2000l) . Thus the Herschel 
Space Observatory {Herschel) and the Planck Surveyor will play 
a crucial role in accurately determining the Rayleigh- Jeans tail of 
the dust SED. The SWIRE XMM-LSS and LH fields have been 
observed as part of the HerMES programme, which will allow us 
to accurately determine the shape of the Rayleigh- Jeans tail of the 
dust SED and provide an insight into the role of cold dust. 



4.5.1 LinktoSMGs 

In this section we plot our 8 2 70 /jm sources on the luminosity- 
temperature (L-T) diagram o f lYangetaljj2007l) to test the link be- 
tween local and intermediate redshift galaxies to the high-redshift 
SMGs. The luminosities are estimated by integrating Eq[T] be- 
tween 40 and 1000 ^.m which we define as Lfir. The two 
galaxy populations appear to lie in different parts of the L-T plane 
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Figure 12. Left: Dust temperature distribution estimated using Eq[T] Right: Dust mass distribution estimated from Eq|2] 



such that the SMGs appear to have colder dust temperatures and 
higher dust masses which is seen as evidence for evolution in 
the dust properties of galaxies. In Fig[T3]we plot the positions of 
our 70 sources alon gside SMGs from Chapman et al. (200^; 
ICoppin et alj ( |2008|) andlKoyacs et alj l2006h as well as moderate- 
z UL IRGs from Yang et al.l ( 120070 ■ more normal SLUGS sources 
from iDunne et alj j200oir local ULIRGs from IClements et"al] 
fcOlOb*). We also plot the 70 ^.m sample of S09 selected at z < 1.2, 
S70 > 9mJy and Sieo > 50mJy by estimating Td and Lfir using 
our method. The main difference between the two 70 /im samples 
is that our work is restricted to sources with spectroscopic redshifts 
while some of the sources in the S09 sample have photometric red- 
shift estimates. Furthermore, only 62% of the S09 sample has > 5(t 
detection at 160 /im. Thus the error estimates for the S09 sample in 
Td and Lfir, are expected to be larger. 

We find that our 70 sources (red stars) in Fig[T3]lie sys- 
tematically below the plane of the SLUGS sources (plus sign) and 
local and moderate-z ULIRGs (crosses and filled circles) in the L- 
T diagram. Our 70 pm sources span a range of temperatures ~ 24 
- 60 K with ~ 30% having Td < 30 K, whereas only ~ 6% of 
the SLUGS sample and none of the moderate and local ULIRGs 
have Id ^ 30 K. This implies that our 70 sources are notice- 
ably cooler than what has been observed previously in the local 
to intermediate IR luminous sources. The 70 ^m sample of S09 is 
characterised by dust temperature Td = 34.3 ± 7.45 K, where ~ 
3 1 % have Td < 30K, which implies that the dust properties of the 
two 70 /im samples are similar. Thus it appears that the 70 /im pop- 
ulation is beginning to fill the gap between the cold high-z SMGs 
and the warmer local to moderate-z IR galaxies. 

The detection of cooler dust can be interpreted as the 70 
sources r epresenting a missing lin k between the two population of 
galaxies jSymeonidis et alj|2009r ) or as a result of selection effect 
jYang et al .120071) . The studies of the local to moderate-z IR galax- 
ies are based on IRAS 60 /xm derived samples and therefore these 
studies may have been biased to warmer dust temperatures whereas 
the high-z SMGs are may be biased towards cold dust tempera- 
tures. Thus using observations at the longer 70 band, we are 
most likely observing galaxies that have cooler dust temperatures 
when compared to the IRAS selected sample and therefore bridging 
the separation between the high-z SMGs and local sources. 

Results from a Pl anck study of nearby galaxies 
jPlanck Collaboration! I20T1I) has found evidence for cold dust 
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Figure 13. Temp erature-luminos ity diagram comparing ULIRGs (crosses: 
|Yang etalj |2007!; filled circles : IClements et alJl2010M . SMGs (squares: 
Chapman et al. 2005; triangles: lKovdcser^l2006l) , normal SLUGS galax- 
ies I D unne et alJ200d) and S09 70 ^tm sample based on single temperature 
dust SED fits, with j3 = 1.5. LpjR i s derived by integrating Eq[T]between 
40 - 1000 fjLm. Diagram adapted from lClements etalj<2010bl) . Typical error 
bars are shown at the bottom right. 



at T < 20K with a variable emissivity and therefore reducing 
the distinction between the two populations. The Planck sources 
overlap with the SLUGS galaxies but extend to colder temperatures 
with temperatures ra nging from 15 - 50K. Studies carried out by 
lAmbIardetaD(l2010l) using Herschel observations have also begun 
to fill in the gap between lo cal IRAS galaxies and high-z SMGs. 
IPlanck Collaboration! ( 1201 ih conclude that cold dust is significant 
and remains an unexplored component in many nearby galaxies. 
As mentioned in the previous section we will be able to further 
understand the nature of the 70 /im sources using observations 
from Herschel, which will provide data points across wide range 
of wavelengths from 70 - 500 /im. 
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5 SUMMARY AND CONCLUSIONS 

We present spectroscopic follow-up observations of 70 pim selected 
sources in the SWIRE XMM-LSS and LH-ROSAT regions using 
the multi-object fibre spectrograph AF2AVYFF0S on WHT. We 
measured new spectroscopic redshifts for 293 70 /im and 35 24 
^m sources. The redshift distribution for the 70 /im sources peaks 
at 2; ~ 0.3 and has a high redshift tail out to 2; ~ 3.5. The majority 
of the 24 fj,m selected QSOs are at typically z > 1. The spectro- 
scopic redshifts of our sample were compared with the photometric 
redshift estimates which show a good agreement for galaxies while 
some QSOs have poorer agreement which we attribute to optical 
variability or photometric redshift aliasing. 

We carry out emission line diagnostic for 91 70 fim sources 
with [OIII], H/3, [Nil], Ha and [SII] emission lines present to clas- 
sify the galaxies into star-forming, Seyfert, composite and LINERs. 
We find in our sample 34 star-forming, 30 composite, 1 Seyfert, 5 
LINERs and 21 ambiguous galaxies and 13 QSOs identified from 
their broad emission line spectra. 

We then modelled the SEDs for each source from optica l to the 
FIR using the method described in Rowan-Robinson et al. II200I 
I2OO8I) and SED templates from lRowan-Robinson et ai]( l2008l) .Our 
major finding is that the IR emission in the 70 /im sources are pow- 
ered by star formation, with the AGN du st torus component mak- 
ing a small contribution for non QSOs. iSvmeonidis et al.l ( 1201 Ol) 
examine the AGN content of a sample of IR luminous 70 /im se- 
lected galaxies by modelling the panchromatic SEDs from X-ray 
to IR and find for most of their sources, the AGN contributes less 
than 10% to the IR budget. Moreover, they find all sources in their 
sample are primarily powered by star-formation . Our results are 
also consistent with the findings of iTrichas et al.l ( |2009|) who find 
their X-ray selected sample with 70 detection are strongly star- 
forming and require a starburst SED to fit the 70 /im and 160 nm 
(if available) photometry. For 9 QSOs and the 1 Seyfert galaxy the 
AGN dust torus dominates in the NIR and MIR regimes produc- 
ing from 40 to 100% of the total IR luminosity. Of the 83 70 fim 
sources in our sample with 160 /im detections, 49 require the cir- 
rus component which implies the presence of large amounts of cold 
dust at Td ~ 25 K. By integrating the SEDs from 8 - 1000 fim, we 
find our 70 fim sample is dominated by starbursts and LIRGs. We 
also find that the AGN fraction increases with Lir with almost all 
(7/10) HLIRGs being optically identified as QSOs. 

Finally, we fit single temperature modified blackbodies to 82 
70 fj,m sources with 160 fim fluxes at 2 < 1.2 with /3 fixed at 1.5. 
Our sample has a mean and standard deviation dust temperatures, 
Td = 32.7 ± 6.7 K and dust mass range, M^ = 10*^ - 10^ Mq and 
even increasing /3 to 2 results in slightly lower temperatures, T^ 
= 30 ± 6 K. In general we find that as /3 increases T^ decreases. 
Examining our sources in the L-T diagram shows that the 70 ^m 
sources with 160 /xm detection have dust temperatures that are sys- 
tematically lower than the local IR luminous galaxies. Herschel and 
Planck will play a vital role in characterising the IR galaxy popu- 
lation, as well as breaking the degeneracy between /3 and Td by 
providing data points across range of FIR/submm wavelengths. 
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